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Abstract 


Two potential coefficient fields that are complete to degree and order 360 
have been computed. One field (OSU86E) excludes geophysically predicted 
anomalies while the other (OSU86F) includes such anomalies. These fields were 
computed using a set of 30' mean gravity anomalies derived from satellite 
altimetry in the ocean areas and from land measurements in North America, 
Europe, Australia, Japan and a few other areas. Where no 30* data existed, 
I'xl* mean anomaly estimates were used if available. No rigorous combination 
of satellite and terrestrial data was carried out. Instead we took advantage of 
the adjusted anomalies and potential coefficients from a rigerous combination 
of the GEML2’ potential coefficient set and l*xl* mean gravity anomalies. 

The two new fields were computed using a quadrature procedure with 
de-smoothing factors suggested by Colombo. The spectra of the new fields 
agree well with the spectra of the fields with l*xl* data out to degree 180. 
Above degree 180 the new fields have more power. The fields have been 
tested through comparison of Doppler station geoid undulations with 
undulations from various geopotential models. The agreement between the two 
types of undulations is approximately *1.6 m. The use of a 360 field over a 

180 field does not significantly improve the compariosn. Instead it allows the 
comparison to be done at some stations where high frequency effects are 
important. In addition maps made in areas of high frequency information 
(such as trench areas) clearly reveal the signal in the new fields from degree 

181 to 360. 
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1. Introduction 


This report describes our first attempt at carrying out a spherical 
harmonic expansion of the earth’s gravity field to degree and order 360. The 
procedure followed is built on an adjustment of satellite derived potential 
coefficients and l*xl*.mean gravity anomalies that is described by Rapp and 
Cruz (1986). This combination solution yielded a set of potential coefficients 
corresponding to those values with a priori informationi and a set of adjusted 
l*xl* anomalies. These anomalies were then used in a least squares collocation 
estimation procedure to determine a set of potential coefficient to degree 250. 
Examination of the spectrum of this field led to the suggestion that smoothing 
was taking place at the higher degrees due to the averaging in the l*xl* 
cells. 

In order to estimate more reliable high degree fields we decided to use 
anomalies in smaller block sizes — specifically 30’x30’ anomalies. A rigorous 
repeat of the prior adjustment was not possible so that an alternate procedure 
was developed to obtain an expansion to degree 360. The following sections 
describe the methods used and the results obtained. 


2. The Previous Combination Solution 


In this section we briefly review the procedures that were followed by 
Rapp and Cruz (1986). The starting equation is the expression for the 
gravitational potential written in the following form: 


V(r, + ,X) = ~ [l ■'■iijlf) + S£mSinmX)Piin(sin*) 


( 1 ) 


where: 


r.^.X 


kM 


a 


Sim 

Pim 


geocentric coordinates 
geocentric gravitational constant 
equatorial radius of the reference ellipsoid 
fully normalized potential coefficients 

fully normalized Legendre function of degree i and order m. 


The disturbing potential can be written as: 
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(3) 


The gravity anomaly, in the Molodensky surface free-air anomaly sense is 
given as: 


Agp = gp “ ro 


(4) 


where P is at the surface and Q is on the telluroid. The boundary condition 
that relates Ag and T is as follows: 


Ag = - 


il + i ±r 

3h 7 3h 


T 


(5) 


where h is the plumbline direction. Neglecting deflections of the vertical 
equation (2) can be substituted into (5) to obtain: 




i i 
£ 

m=o 




(6) 


where (i=h, 7 ) are ellipsoidal corrections. If we assume that we have the 

anomalies (Ag) given on the surface of the reference ellipsoid equation (6) can 
be inverted to solve for potential coefficients (Cruz, 1986). We have: 


= Cii°+ Cl^+ Cii"+ 


(7) 


In (7) we have: 


Cii°= ^ ^ I I AgEYfm(*,X)d<r 
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( 8 ) 
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The other terms in (7) are ellipsoidal correction terms defined in Rapp and 
Cruz (ibid, pp. 16-17). 

The gravity anomaly on the ellipsoid must be determined by downward 
continuing the surface anomaly to the ellipsoid and adding an atmospheric 
correction (dg^)* We have used the following: 


= Ags + <5gA - 


h - 


A 

2! 3r^ 


h^ 


(9) 


This downward continuation is an approximation. For our work the derivatives 
were calculated from an expansion to degree 300 which was developed as part 
of a previous solution (Rapp, 1981). 

The combination solution was carried out through the adjustment of the 
potential coefficients derived from satellite analysis and the coefficients 
implied by gravity data. Specifically the mathematical model, F, was taken as 


F — l*x^ "" a — 0 (10) 

c 

where Lxa. is the adjusted potential coefficient set derived from satellite data 

and Lx ^ potential coefficient set implied by the adjusted gravity 

c 

anomalies given in cells on the earth. The linearization of (10) is written; 


BiVi + BxVx + W = 0 


( 11 ) 


where V| are corrections to the approximate anomalies, and Vx are corrections 
to the a priori known potential coefficients. The W vector is: 


W = Lxo - Lx^c (12) 

where LxO is the approximate potential coefficient set and Lxq® is the 

potential coefficient set implied by "observed” anomalies. The model is 
formulated taking into account ellipsoidal corrections and downward 

continuation effects. We start this by computing the values from the 

given potential coefficient using (7). We have: 


= LxO = Cin - 


( 13 ) 


where Acfi® are all but the first terms on the right side of (7). Using 
equation (8) the values of Lx^c s^e computed. 
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The adjustment is now carried out to obtain and Vi. Specifically: 


Vx = - ((Bi Pi^ bJ)-» + Px)-‘ (it PV bI)W (14) 

yt = PJ‘ b] Px Vx (15) 


Pi is the weight matrix of the observations and Px is the weight matrix of the 
a priori potential coefficients. The elements in Bi are given by the 
coefficients of Agg in equation (8). The adjusted potential coefficients, 
compatible with satellite derived estimates, would be: 


Cim,a = Lxo + Vx + ACii° 


(16) 


The adjusted anomalies on the ellipsoid would be: 


Lia = L^o + Vi 


(17) 


while the adjusted anomalies at the surface of the earth would be found from 
equation (9). 

In our combination solution we used the GEML2' potential coefficients that 
are very similar to the GEML2 coefficients described by Lerch et al. (1982). 
In addition we used a few additional coefficients to degree and order 30. The 
total number of coefficients estimated was 582. 

The I'xl* gravity anomaly data was based on the merger of terrestrial 
anomaly estimates with altimeter derived gravity anomalies. In this merger we 
considered two forms. One form excluded most geophysically predicted 
anomalies while the second form included such anomalies. In the first case we 
had a total of 50,562 l*xl* anomalies while in the second case there were 
56,109 values. Figures 3 and 4 of Rapp and Cruz (ibid) show the location of 
these anomalies. In order to fill in the "empty” areas we computed the 
anomalies from the a priori potential coefficient set. 

The higher degree potential coefficient fields were made by applying an 
optimal estimation procedure to the adjusted anomalies on the ellipsoid. This 
would correspond to a least squares collocation estimation using appropriate 
covariance functions, anomaly degree variance models, and anomaly accuracy 
estimates. The process of the optimal estimation as developed by Colombo 
(1981) with implementation for l*xl* data described by Hajela (1984) requires, 
for a complete solution, extensive computations. The solution we estimated was 
to degree 250. After the coefficients are estimated it is necessary to add the 

corrections introduced in equation (13). 
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It is possible to approximate the optimal estimate through the use of a 
modified orthogonality procedure. If the anomalies were given on the 

complete surface of the ellipsoid we could use ( 8 ) to compute 0 ?^° • Since 
data are actually given in discrete cells an implementation of ( 8 ), to 
approximate an optimum estimation, was suggested by Colombo (1981). If the 
data is given in blocks of 9* let N = 180*/9*. Then we have: 

4ny(i-l)q, % I I (18) 

where: 7 = kM/a^; 

Agjj = the mean anomaly corresponding to Agg; 

qg = fit; 0 ^ t ^ N/3 
qi = fit; N/3 < A i N 

qi = 1 ■« > N (19) 


fijt is the Pellinen smoothing operator that depends on a circular cap radius 
which in turn depends on the size of the mean anomady cell (Rapp and Cruz 
(ibid, p. 22 )). 

Using the adjusted l*xl* mean anomalies of our combination solution we 
computed a set of potential coefficients from equation (18) which we have 
compared to the values from the optimal estimation solution. Table 1 (taken 
from Rapp and Cruz (ibid. Table 11) shows the differences in terms of 
undulations, anomeJies, and percentage differences. 


Table 1. Comparison of the Potential Coefficient Solution from Equation (18) 
with the Optimal Estimate Using l*xl* Mean Anomalies. 


t 

1 5N (cm) 

Ifig(mgals) 1 

PCX) 

2 

0.0 

0.00 

0,0 

10 

0.0 

0.00 

0.0 

20 

0.0 

0.00 

0.0 

30 

0.2 

0.01 

0.5 

50 

0.2 

0.02 

1.1 

75 

0.8 

0.09 

6.2 

100 

1.0 

0.15 

10.6 

120 

1.0 

0.18 

13.2 

150 

0.6 

0.14 

13.0 

180 

0.8 

0.22 

22.3 

to 180 

8.9 

1.53 

7.6 
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We see from this table that there is good agreement between the two solutions. 
The percentage differences are negligible at the lower degrees increasing to 
22% at degree 180. The results from this test will be used later to justify the 
use of equation (18) with 30’ mean anomalies instead of using the complicated 
optimal estimation results. 

This section has been written to describe our potential coefficient 
solution using l*xl* anomaly data and the GEML2’ potential coefficients. The 
results of these computations were two potential coefficient sets computed 
using an optimal estimation procedure with a uniform data noise of one mgal 
for the anomalies. The solution that excludes geophysically predicted 

anomalies was designated OSU86C; the solution that included such anomalies 
was called OSU86D. Both solutions were complete to degree and order 250. 


3. The 30*x30* Mean Anomaly Data Set 

In order to extend the OSU86C/D solutions to a higher degree we first 
needed to put together a set of smaller anomaly cells in order to get the 
higher frequencies into our solution. The cell size chosen was 30’x30’ 
although other sizes could be chosen. However this size is convenient to 
work with as four 30’x30’ cells fill a I'xl* cell. 

The 30’x30’ mean anomalies for the ocean areas had been computed by 
Rapp (1986) from Geos-3/Seasat altimeter data. Specifically computed were 
point anomalies on a 1/8* grid which were meaned to compute both 30’ and 1* 
mean anomalies (and sea surface heights). The total number of predicted 30' 
anomalies was 149, 670 although about 1000 were considered unreliable because 
they were on land or in areas in which the altimeter data was sparse. 

In order to obtain data from terrestrial estimates it was necessary to 
collect gravity anomaly data as point values or small (lO’xlO’, 6’xlO’) mean 
values. Such values were then formed into 30’x30’ values with individual 
accuracy estimates where possible. The development of this 30’ terrestrial 
anomaly data base is described by Despotakis (1986). Figure 1 shows the 
location of the 30’ mean anomalies in the August 1986 field. 

We next created a merger of the altimeter derived anomalies and the 
terrestrial values. In the ocean areas we used the altimeter derived anomaly 
unless its standard deviation was greater than 30 mgals or if the anomaly was 
in a l*xl* block whose elevation was greater than or equal to zero. For the 
remaining areas the 30’ terrestrial values were used if they were available. 
The atmospheric correction (6g/^) was added to the terrestrial anomalies (see 
equation 9) and a gravity formula correction of -0.6 mgal was added to 
convert to our reference system of constants. If neither an altimeter derived 
anomaly nor terrestrial anomaly was available the adjusted I'xl* anomaly from 
the combination solution was used as fill-in values. 
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Figure 1. Location of 31450 30’x30 Mean Free Air Anomalies Estimated From Terrestrial 
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with at least a single 30 * value from the terrestrial or 








In order to take advantage of the least squares adjustment that led to the 
OSU86C/D solutions we tooH the adjusted l*xl* values of these solutions and 
forced the mean of the four 30' values within the blocks to be the same by 
applying a bias term to each of the non-fill-in 30' values. Specifically we 
write: 


B = 4gE(l*.A) 


SAg^OO'.UA) 

n 


( 20 ) 


In (20) the A stands for adjusted and UA stands for unadjusted and n is the 
number (usually 4) of non-fill-in 30' values in the 1* blocks. The B term 
would consist of three effects: the downward continuation of the surface 

anomalies to the ellipsoid (see equation 9); the inconsistencies of the 30' and 
1* values; and differences caused by the adjustment process. The B value for 
each 1* block was then added to the non-fill-in 30' value to obtain a value 
that we consider to be on the ellipsoid (Agg) and is consistent with the 
adjusted anomalies of our earlier solutions. 

In summary the merger process used 139,946 values derived from altimeter 
data; 21,739 values from the terrestrial data, and 97,515 "fiU-in" values from 
the adjusted l*xl* data. The total number of 30' anomalies is 259,200. Figure 
2 shows the location of 40700 l*xl* blocks in which at least one 30' value from 
either the altimeter derived or terrestrial data set was used. 


4. The Estimation of the Potential Coefficients 


If we were to start the combination process from our 30' data we would 
need to evaluate equations (14) and (15). This process is feasible with a 
vector processing machine. However we did not feel that the effort for such 
a solution was worthwhile at this time. Since the actual adjustment is done 
with the coefficients of the satellite solution, which was complete only to 
degree and order 20, we felt that whether this adjustment was made with 1* 
data or 30' data would be irrelevant. We consequently adopted the ad justed 
coefficients of the OSU86C/D solutions. For the other coefficients to degree 
360 we used equation (18) through the HARMIN program of Colombo (1981) 
where the Ag j j values were taken as the bias corrected 30' values. This 
computation required an evaluation of the integrals of the fully normalized 
associated Legendre functions up to degree and order 360 at 30' intervals. 
This was done using program F428AV1 which implements the Paul (1978) 
subroutine with changes suggested by Gleason (1985). 

The HARMIN program (OSU program F419B) was used with the qi factors 
given by equation (19). Two solutions to degree 360 were carried out. The 
OSU86E solution basically excluded geophysically predicted anomalies while the 
OSU86F solution included such anomalies. After the potential coefficients 

from equation (18) were fovmd, the ellipsoidal corrections given by AC^i® of 
equation (13) were applied. These corrected coefficients were then merged 
with the adjusted and corrected coefficients of the OSU86C or D solution to 
form the final potential coefficient sets. 
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5. Potential Coefficient Accuracy 


One disadvantage of using the HARMIN approach to the potential 
coefficient estimation is the lack of estimated coefficient accuracies. In fact, 
the use of the more rigorous optimal estimation procedure does not ensure 
realistic accuracies because of the smoothing of the field that takes place at 
the higher degrees when realistic anomaly standard deviations are used. To 
avoid this problem Rapp and Cruz (ibid) combined the sampling error (from 
the optimal estimation procedure) with a propagated noise assuming 
uncorrelated anomaly errors of 10 mgal. This procedure led to a noise 
magnitude that was equal to the signal near degree 175 (see Figure 3). Since 
we do not have a formal sampling error estimate for the 30' solution, and 
since the assumption of uncorrelated 30' anomaly errors is even more 
unrealistic than the 1* case we have not computed the potential coefficient 
errors for the OSU86E and ’ F solutions. If there are requirements that 
demand coefficient accuracy estimates for the 86E/F fields, we suggest the use 
of the 86C/D standard derivation to degree 175 followed by a 100% uncertainty 
in the coefficients above that degree. 


6. Anomaly Degree Variances 

We first define the spectrum or power at degree i- as follows: 


(Tl“ O 


( 21 ) 


The anomaly degree variances, are given as; 


Ci = (22) 


where y - kM/a*. The anomaly degree variances are formally given on a 
sphere of radius a. Figure 3 shows the anomaly degree variances for the 
OSU86C and the OSU86E solution. Also shown is the accuracy estimate of the 
OSU86C solution based on the optimal estimation procedure with an anomaly 
standard deviation of one mgal and a propagated standard deviation of 10 
mgal for the l*xl* anomalies (Rapp and Cruz, ibid, equation (6.11)). We see 
that the power is almost identical in the two solutions out to degree 200 
although OSU86C solution shows slightly less power between degrees 120 and 
160. Beyond degree 204, out to degree 250, the 86C solution has power that 
falls off much more rapidly than the solution (OSU86E) with the 30' anomalies. 
This seems to reflect the loss of information, not recovered by the optimal 
estimation procedure, in the 1* averaging process. This might imply that the 
rough rule of thumb (18O*/0*, 9* = block size) for how high a degree should a 
spherical harmonic expansion be taken may not be unreasonable. 
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Figure 4. Anomaly Degree Variances for OSU86D (1*, with geophysical anomalies) and OSU86P (30», with 
geophysical anomalies); Error Anomaly Degree Variances for OSU86C/D. 



Figure 3 also shows the anomaly degree variance implied by the Kaula rule 
for the decay of potential coefficients (i.e. the rule). The anomaly 

degree variance implied by this rule has been computed on a sphere of radius 
a assuming the original rule was based on data referred to the mean surface 
of the earth. Specifically plotted was: 


ci(a) = 


192 

i+1.5 



n+2 


mgal* 


(23) 


where Rg = 6371 km and a = 6378 km. The rule has too much power between 
degrees 10 to 50 and too little power between degrees 80 to 200. Beyond 200 
the estimates from the model and from OSU86E are within 50% of each other. 

Figure 4 shows the anomaly degree variances for the OSU86D (1*) and the 
OSU86F (30’) solutions. Both these solutions used the geophysically predicted 
anomalies. Also shown is the accuracy of the OSU86C/D solutions. The 
comments made with respect to the anomaly degree variances shown in Figure 
3 are also applicable here. 

Specific values of the anomaly degree variances, at selected degrees, for 
the four solutions discussed here are given in Table 2. 


Table 2. Anomaly Degree Variances. (Units are mgaP). 



Solution 

i 

0SU86C 

0SB86E 

0SU86D 

0SU86F 

50 

2.66 

2.57 

2.94 

2.83 

100 

2.23 

2.13 

2.35 

2.23 

150 

1.36 

1.28 

1.48 

1.33 

180 

0.96 

0.97 

1.02 

1.01 

200 

0.57 

0.60 

0.56 

0.61 

250 

0.14 

0.41 

0.16 

0.43 

300 


0.22 


0.22 

350 


0.18 


0.18 


From this table we see that there is no substantial difference, at the higher 
degrees, in the power of the OSU86E and F solutions. 


7. Doppler Undulation Comparisons 

As discussed by Rapp and Cruz (ibid, section 8.3) the accuracy of a 
potential coefficient model may be judged by comparing the geoid undulation 
derived from the model and the vsdue implied by the ellipsoidal and 
orthometric heights at a station. This can only be done after the station 
coordinates have been converted to a geocentric, true scale system, and the 
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best estimate of the equatiorial radius is used. 

We carried out these computations with the new models for a global 
Doppler station set, and data sets in North America, Australia and in Europe. 
In these tests an equatorial radius of 6378136 m was used. the mean 

difference (Doppler minus model) and the standard deviation of the difference 
is given in Table 3 for Europe, and Table 4 for Australia for a fixed number 
of stations. 


Table 3. Comparison of Doppler Derived Undulations in Europe With Values 
From Models. 173 Stations Used. 



Table 4. Comparison of Doppler Derived Undulations in Australia With Values 
From Models. 114 Stations Used. 



From these two tables we see that in Europe, of the degree 180 solutions the 
OSU86E is best, while in Australia GPM2 is the best. The degree 360 fields in 
Europe show a slight improvement over their 180 counterparts. In Australia 
the 360 fields are slightly better than the 180 fields with the exception of 
GPM2 which fits the best for this area. 

For North America and the global data set we carried out the comparisons 
rejecting stations where the residual was greater than 4 m. Therefore in 
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judging the results of these comparisons we consider the standard deviation 
and the number of stations accepted for the comparison. These results are 
shown in Table 5, for North America, and in Table 6, globally. 

Table 5. Comparison of Doppler Derived Undulations in North America With 
Values from Models 


Model 

i(max) 

Mean 

Difference 

Std Dev 
Difference 

Num of 
Stations 

0SU81 

180 

0.22 

*1.72 m 

683 

GPM2 

180 

0.18 

1.58 

695 

0SU86C 

180 

0.20 

1.54 

687 

0SU86D 

180 

0.26 

1.55 

687 

0SU86E 

180 

0.21 

1.52 

689 

0SU86E 

360 

0.26 

1.45 

695 

0SU86F 

180 

0.28 

1.52 

688 

0SU86F 

360 

0.31 

1.47 

696 


Table 6. Comparison of Doppler Derived Undulations Globally With Values From 
Models. 


Model 

i (max) 

Mean 

Difference 

Std Dev 
Difference 

Niam of 
Stations 

0SU81 

180 

0. 18 ffl 

*1,76 m 

1721 

GPM2 

180 

0.13 

1.64 

1735 

0SU86C 

180 

0.17 

1.66 

1741 

0SU86D 

180 

0.21 

1 . 66 

1743 

0SU86E 

180 

0.16 

1.67 

1752 

0SU86E 

360 

0.17 

1.66 

1777 

0SU86F 

180 

0.21 

1.66 

1754 

0SU86F 

360 

0.21 

1.66 

1780 


From Table 5 we see, of the 180 solutions, both the OSU86E/F solutions 
give essentially the same results. Improved results are found with the 360 
solutions with a small decrease in the standard deviation of fit and a small 
increase in the number of stations accepted. 

From Table 6 we see all the solutions except the OSU81 solution give about 
the same standard deviation of fit. Of the 180 solutions, the OSU86F solution 
accepts the most stations. Of the two 360 solutions tested the OSU86F solution 
accepts the most stations. 

In judging these results we must recall that the accuracy of the Doppler 
derived undulations may not be sufficient to distinguish the accuracy of 
various models at or below the *1.5 m level. More accurate results might be 
obtained using laser station coordinates. Unfortunately such stations are 
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much fewer than the Doppler stations and they are located at difficult 
locations (e.g. mountain tops). 

Rapp and Cruz (ibid, Section 8.4) discussed the undulation residual 
correlation with topography. Such correlation continues for these 360 fields. 
For example for the OSU86F field to degree 360, with the global station set, 
the slope is (0.41*. 14) m/km with 1780 stations accepted. The value for the 
OSU86D field to degree 180 is (0.41*. 12) m/km with 1741 stations accepted. 


8. High Degree Geoid Maps 

It would be possible, but not really meaningful, to compute a global geoid 
undulation map using the two new potential coefficient models to degree 360. 
Such maps would not be meaningful because the high degree contributions to 
the undulation, for the most part, are quite small. To first estimate the 

magnitude of the high degree terms we computed, for selected degrees, the 
undulation contribution from degrees i+1 to degree -^max* Values for selected 
degrees are given in Table 7. 

Table 7. Geoid Undulation Magnitudes for Selected Degree Ranges in the 
OSU86F Model 


Degree Range 

Undulation Contribution 

180 to 250 

*20.2 cm 

181 to 360 

23.2 cm 

251 to 300 

9.3 cm 


The values given in Table 7 represent global averages. To see specific 
estimates of the contribution made by the higher degree terms, we prepared 
two maps showing the geoid undulations computed from the coefficients from 
degree 181 to 360. Figure 5 shows this map for a region in the United States 
covering a portion of the Rocky Mountains. Figure 6 shows these effects for 
a region over the Tonga Trench. 

In the Rockies test we can see effects that reach 1.8 m. In the Tonga 
Trench area the high frequency signal over the trench is clearly visible in 
the lineated pattern. The largest values are on the order of 2 meters. The 
latter case shows how the high degree fields can play a role where high 
frequency information is present in the gravity field. Note from Figure 6 
that, outside the trench the contributions of the higher frequencies are quite 
small, approaching, it is speculated, the *23 cm computed for Table 7. 

Forsberg (1986) has computed the geoid undulation information above 
degree 180 based on power spectral analysis of gravity anomalies in Norway, 
Sweden, Finland and Denmark. Averaging over 37 2*x4* areas the information 
above degree 180 is *32 cm with uncorrected anomaly data and *28 cm with 
terrain corrected anomaly data. These values compare well with the *23 cm 
value given in Table 7 when the summation is taken just to degree 360. 
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Figure 6. Geoid Undulations in the Tonga Trench Area from Degrees 181 to 
360 of the OSU86E Solution. The contour interval is 0.5 m. The 
data grid was 0.*25. 












9. Coefficient Comparisons 


It is of interest to compare the potential coefficients of the C/D solution 
with coefficients of the E/F solution up to degree 180. The C vs E comparison 
is given in Table 8 while the D/F comparison is given in Table 9. 

Table 8. A Comparison of the OSU86C and E Models at Selected Degrees 


i 

5N(cnj) 

5<?(aigals’) 

Percent Diff. 

50 

0.7 

0.05 

3.3 

75 

0.9 

0.11 

6.9 

100 

1.2 

0.18 

11.9 

120 

1.4 

0.26 

17.2 

150 

1.7 

0.38 

33.0 

180 

1.8 

0.49 

50.2 

2 to 180 

15.7 

3.18 

14.6 


Table 9. A Comparison of the OSU86D and F Models at Selected Degrees 


_£ Syi(an) ^g(mgals) Percent Diff. 

3.5 

6.9 

11.7 

17.9 

32.8 

48.9 
14.6 


50 

0.8 

0.06 

75 

1.0 

0.11 

100 

1.2 

0.18 

120 

1.5 

0.27 

150 

1.7 

0.40 

180 

1.8 

0.49 

to 180 

16.1 

3.25 


The main thing that we see from these two tables is that the differences 
between the 1* and 30* data solutions are small up to degree 180. The 
differences are a function of degree, with the largest differences occuring at 
the higher degrees. This gives us some confidence in the 1* solutions but it 
also shows us the coefficients at the higher degrees may change by 50% from 
a 1* solution to a 30’ solution. It would not be unreasonable to expect the 
coefficients at the high degrees of the E/F fields to change by this amount if 
0.*25 anomalies were used instead of 0.*5 values. 

As a last comparison we show in Figure 6 the undulation differences 
between the OSU86F and the OSU81 solutions up to degree 180. The contour 
map has been created from data given on a 2*x2* grid so that some high 
frequency differences may be missing. The maximum difference is 9.7 m which 
occurs in the southwest part of Africa. Most of the large undulation changes 
are due to substantial anomaly differences in the solutions. 
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Figure 7. Geoid Undulation Differences (OSU86P - OSU81) Contour Interval = 1 m; 2*x2* Grid. 




10. Conclusions 


This report has described the estimation of an expansion of the earth’s 
gravitational potential to degree 360. The method used, did not rigorously 
combine, satellite and terrestrial gravity data. Instead, we built on solutions 
that were formed using l*xl* anomalies. We used a 30’ data base formed from 
the merger of the anomalies derived from satellite altimetry in the oceans and 
terrestrial data for some land areas. Many land areas have no coverage or 
coverage only in terms of l*xl* anomaly estimates. In such areas the field 
will not represent the high frequency variations that may actually exist in the 
area. 


An important part of the anomaly reduction process is the downward 
continuation of the surface anomalies to the ellipsoid. In dealing with the 1* 
data we used a two term Taylor series involving a preliminary high degree 
field. We did not compute the downward continuation effects for the 30’ cells. 
Instead we effectively used for each 30’ cell the value used in the 1* 
reductions. Future solutions should examine improved techniques to carry out 
this downward continuation process. 

We did not compute accuracy estimates for the coefficients. Such 
estimates would require a number of assumptions that would be unrealistic in 
practice. More work needs to be done in this area. 

The next generation of the 360 field could involve the formal adjustment 
of the 30’ mean anomalies and satellite data. More complete downward 
continuation procedures could be attempted. Accuracy evaluation should be of 
high importance. 

Although there are a number of improvements that could be made for the 
next solution, the current solution appears to represent the given data and 
external data (e.g. Doppler undulations) quite well. We must remember, that 
just because we have a high degree field, it does not mean we have a highly 
accurate high degree field. The latter type of field will only be available 
when more precise data, such as would come from the Geopotential Research 
Mission of NASA (1986) would become available. Until then our solutions to 360 
can be used for a number of purposes including realistic simulations for the 
GRM mission. 
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